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Counterion distribution around an isolated flexible polyelectrolyte in the presence 
of a divalent salt is evaluated using the adsorption model [M. Muthukumar, J. Chem. 
Phys. 120, 9343 (2004)] that considers Bjerrum length, salt concentration, and local 
dielectric heterogeneity as physical variables in the system. Self consistent calcula- 
tions of effective charge and size of polymer show that divalent counterions replace 
condensed monovalent counterions in competitive adsorption. The theory further 
predicts that at modest physical conditions, polymer charge is compensated and 
reversed with increasing divalent salt. Consequently, the polyelectrolyte collapses 
and reswells, respectively. Lower temperatures and higher degrees of dielectric het- 
erogeneity enhance condensation of all species of ions. Complete diagram of states 
for the effective charge calculated as functions of Coulomb strength and salt con- 
centration suggest that (a) overcharging requires a minimum Coulomb strenth, and 
(b) progressively higher presence of salt recharges the polymer due to either elec- 
trostatic screening (low Coulomb strength) or negative coion condensation (high 
Coulomb strength) . A simple theory of ion-bridging is also presented which predicts 
a first-order collapse of polyelectrolytes. The theoretical predictions are in agreement 
with generic results from experiments and simulations. 
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I. INTRODUCTION 



Understanding charged polymers has again become a major focus of the polymer com- 
munity in the last few years. With respect to uncharged systems, there are two additional 
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length scales traditionally considered in theoretical treatments of salty polyelectrolyte so- 
lutions. One is the Bjerrum length lg, which sets the length scale for the strength of the 
Coulomb interaction at a particular temperature T in a specific solvent of dielectric constant 
e, given by 



where e is the electron charge, eo is the vacuum dielectric constant, and kg is the Boltzmann 
constant. The other one is the Debye screening length k" 1 , which sets the length scale for 
screening due to dissociated ions, given by 



where the sum is over all species (z) of mobile ions of valency Z$ and concentration q. In 
addition to the above two, a third length scale can be conceived 1 relating to the dielectric 
heterogeneity, which accounts for the difference in the dielectric constant in the vicinity of 
the chains and in the bulk solvent. In this paper, we consider these three length scales to the- 
oretically determine the effective charge and conformation of a single isolated polyelectrolyte 
chain in the presence of divalent salts. 

The typical non-monotonic dependence of the average conformation of flexible polyelec- 
trolyte chains on temperature is now well-known. At very high temperatures, the chains in 
dilute salt-free solutions are in their athermal states with self-avoiding-walk (SAW) statis- 
tics. Although the chains are fully charged at these conditions, the electrostatic repulsion 
among monomers remains negligible compared to thermal fluctuations. As the tempera- 
ture decreases, electrostatics becomes progressively important and the chains expand due to 
inter-monomer repulsion even beyond the excluded-volume swelling. At even lower temper- 
atures, counterions condense on the chains sufficiently reducing the net polymer charge, and 
consequently the chains contract again. The presence of small-molecular monovalent salts 
has long been known to enhance this condensation effect, and the resulting collapse occurs 
at higher temperatures. Monovalent salts, however, collapse polyelectrolytes at tempera- 
tures that are still typically way below modest temperatures ( e. g., room temperature for 
aqueous solutions). Presence of divalent (or multivalent) salts, however, induces a drastic 
qualitative change in polyelectrolyte behaviors. A modest number of divalent counterions 
in water can effectively neutralize and collapse polyelectrolytes at room temperatures. Fur- 
ther, addition of a higher amount of divalent salt can even reverse the charge on the polymer 
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(the phenomenon known as overcharging or charge inversion or charge reversal) at certain 
physical conditions. 

The condensation of counterions on flexible polyelectrolyte chains has traditionally been 
covered in the Manning model 2 originally designed for infinitesimally thin and infinitely 
long rod-like molecules. However, Manning's argument has been found inadequate^ for flex- 
ible polyelectrolytes, particularly for complex systems with multivalent ions. Flexibility 
allows significant bending of molecules due to charge compensation at lower temperatures, 
allowing substantial changes in the conformational entropy of the polymer. Manning's as- 
sumption that the discrete nature of the charged groups has a secondary effect becomes 
entirely invalid for multivalent ions. It is precisely this discreteness that is found respon- 
sible for complete charge compensation (and resulting contraction of polyelectrolytes) and 
subsequent overcharging at modest temperatures by multivalent salt counterions. This over- 
charging behavior is unexplainable within the Poisson-Boltzman formalism which considers 
a continuum description of the charge density. In order to address precipitation of chains for 
high counterion valence, initial theories 3 considered translational free energy of polyions and 
salt ions alongwith screened Coulomb interaction between charges. With prefixed values of 
the excluded volume exponent v (i.e., prefixed radius of gyration, R g ), the free energy was 
minimized in terms of counterion species, and the correlated multivalent ions were shown 4 to 
induce attraction between monomers (through ion-'bridging') capable of collapsing a chain. 
Redissolution of chains was also observed at higher (multivalent) salt concentrations, but 
it was explained by a reduced bridging force due to electrostatic screening (as opposed to 
overcharging). Later, unscreened Coulomb interaction within condensed ion-pairs was first 
addressed^ without considering the chain entropy, and the theory predicted dependencies of 
the degree of ionization, /, (which is the total effective charge density of polyelectrolytes 
after accounting for the condensed ions) on temperature and salt concentrations to be sim- 
ilar to Manning's argument. A two-state (rod-like and collapsed) model for condensation 
predicted 6 that chain collapse occurs when the total charge of multivalent cations equals to 
that of the ionizable groups of the polymer, implying the condensation of almost all added 
multivalent ions at modest temperatures. The two-state theory 6 -^ treats the collapsed state 
at low temperatures as an amorphous ionic solid similar to simple electrolytes (say, NaCl), 
and, therefore, still ignores the chain entropy and bending-related reorganization of con- 
densed charges at low temperatures. 
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Generic experiments^*^ 1 ^ 1 ^ 1 ^ and simulations^ have shown that added cations with 
higher valence are more effective in compactifying a long DNA molecule (or a polyanion in 
general), and that implies a dominant electrostatic mechanism for polyelectrolyte collapse. 
Within the purview of this concept, the issue of overcharging induced by the addition of 
multivalent salts was investigated theoretically for rodlike DNA— molecules. It was 
further noticecr^ 1 ^ that there is a typical range in salt concentration, in which the short- 
ranged attractions between monomers are effective due to the proximity of the isoelectric 
point (at which the average effective charge of the polymer is zero). The total effective 
DNA charge was found to reverse sign if the multivalent salt concentration was increased 
above this range. This concentration window, within which a flexible chain remains col- 
lapsed with virtually negligible net charge on it, is predicted to be very small in recent 
simulations^, which otherwise observe definite overcharging. Although rapid collapse due 
to charge neutralization for higher salt concentrations or lower temperatures is abundantly 
observed in experiments^ 1 ^^ anc j j n simulations^^ 1 ^ 1 ^, charge reversal by multivalent 
salt is not universally observed^. Very recent experiments, however, have lent support to 
both charge inversion^ 1 ^ and counterion mediated attraction^- in biological polyelectrolytes 
such as DNA. 

An effective two-parameter theory^ considered the adsorption process of counterions 
and formation of mono and di-complexes between negatively charged monomers and diva- 
lent positively charged salt ions. The theory predicted charge neutralization and subse- 
quent charge-reversal of the chain backbone at moderate concentrations of divalent salts. 
The predictions, however, are limited to weak polyelectrolytes without consideration to the 
chain energy and hence chain configurations at various physical conditions. To address the 
chain entropy of flexible polyelectrolytes and its role in counterion distribution, Muthuku- 
mar developed 1 a continuum theory of counterion condensation as an adsorption process. 
Condensation in this argument is facilitated at lower temperatures but, unlike in previous 
theories, is coupled with the configurational free energy of the polymer. The adsorption 
theory^ considers continuous values for the size (R g ) of a single chain treated as a continous 
curve, which provide an appropriate description of entropy of a flexible chain. The theory 
treats the salt-free and salty conditions with monovalent counterions in dilute solutions of 
flexible polyelectrolytes. The parametric analysis of the competition bewteen the transla- 
tional entropy of counterions and the electrostatic energy gain of condensed ions reproduces 
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all classical results including the chain-collapse due to short-ranged dipole interactions at 
low temperatures^ 1 ^. In addition to the length scales Ib (Bjerrum length) and k _1 (De- 
bye length) in the charged system, the adsorption theory uses the concept of a dielectric 
mismatch parameter, 5, which captures the fact that the dielectric constant has much lower 
values near the chain backbone of a polyelectrolyte or protein than in the bulk— M£k solvent. 
In generic polyelectrolyte solutions, 5 is the ratio of the bulk to local dielectric constants, 
and the range in which e assumes its bulk value sets a new length scale. This mismatch in e, 
if substantial, will create higher potential gradients that can electrostatically guide counte- 
rions toward oppositely charged monomers. The theory showed that this may significantly 
increase counterion condensation at modest temperatures leading to a lower effective charge 
and smaller size of the polyelectrolyte. Monovalent counterions, however, was shown not to 
be able to collapse a chain completely at modest temperatures (say, the room temperature 
at which Ib ~ 7 A in water). 

In this paper, we extend Muthukumar's adsorption theory^ to include divalent counteri- 
ons. The basic concept still relies on the competition between the electrostatics of condensed 
ions and the entropy of free ions. The present model is intended to analyze the competitive 
displacement of monovalent counterions by divalent counterions when a salt-free dilute so- 
lution of flexible polyelectrolytes is mixed with a salt solution of divalent ions. As a specific 
example of this situation, we consider a single isolated polyelectrolyte (NaPSS) in a dilute 
solution (water) with monovalent counterions of its own but with no additional monovalent 
salt (NaCl) in general. We monitor the chain conformation and arrangement of condensed 
ions as functions of temperature, degree of dielectric mismatch, and concentration of divalent 
salt (BaCy. In addition to the condensation of monovalent and divalent counterions (Na + 
and Ba ++ , respectively), we also take into account the attachment of negative salt coions 
(CI - ), as previous theories^S and simulations^ indicate substantial presence of negative 
coions near the chain backbone. The configurational free energy of the system depends on 
the extent of adsorption of various species of ions on the chain backbone as well as on the 
size of the polymer, and all contributing factors remain non-trivially coupled. These factors 
are assessed self-consistently with the important assumption of the adsorption theory being 
that the chemical potentials in the adsorbed and free states of the ions are the same. A 
numerical minimization of the free energy with respect to the fraction of condensed ions and 
the chain size determines the equilibrium values of the respective quantities. 
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The key conclusions are the following. In a competitive adsorption process, divalent ions 
displace the condensed monovalent ions at modest temperatures and for reasonable values 
of the dielectric mismatch parameter. At similar physical conditions, increasing divalent 
salt concentration can induce complete neutralization and resulting contraction of the poly- 
electrolyte chain. With further increase of divalent salt in the solution, condensed divalent 
counterions overcharge and reexpand chains for significantly large ranges of physical param- 
eters S and Ib- A substantial fraction of negative salt coions (CT~) also condenses on the 
monomer- divalent ion-pair which reduces the degree of overcharging but does not eliminate 
it. However, for modest temperatures in a generic solvent (water at room temperature), there 
is always a critical degree of dielectric heterogeneity below which neither chain collapse nor 
charge reversal occurs irrespective of the amount of divalent salt in solution. Further, the 
dependencies of the regularization of charge and size of the polymer on salt concentration, 
temperature, and the dielectric mismatch parameter are relatively robust functions of Ib, 
and 5, not depending sensitively on the microscopic details of the charge complexes. 
For a moderate presence of salt, both the degree of ionization and chain size are typically 
smaller for divalent salts than for monovalent salts. A higher dielectric mismatch and a 
lower temperature enhance condensation of all types of counterions and coions resulting in 
the achievement of isoelectric point at lower divalent salt concentrations. A typical state 
diagram for polyelectrolyte charge is predicted (Fig. [1]) in which below a critical Coulomb 
strength (proportional to 1/T, 1/e or 8) there is no overcharging with increasing divalent salt. 
Above this critical strength, electrostatics is strong enough to condense sufficient number of 
counterions inducing overcharging. For higher salt concentrations, we expect recharging of 
the polyelectrolyte due to either screening of electrostatics^ (low Coulomb strength) or neg- 
ative coion (Cl~) condensation (high Coulomb strength). The dashed line indicates that in 
this regime of salt concentration the theory only predicts qualitative results. If ion-bridging 
is present, a simple theory based on this model predicts first order chain collapse. Our the- 
ory shows that overcharging is an outcome of both correlation-induced adsorption related to 
the discreteness of divalent cations and dielectric heterogeneity related to the local chemical 
structure of polyelectrolytes. 

The rest of the article is organized as follows: in Sec. II, we develop the theory and 
in Sec. Ill, we present the results discussing competitive adsorption (III. A), chain collapse 
(III.B), overcharging (III.C), coion condensation (III.D), free energy (III.E), state diagrams 
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for various regimes of polymer charge (III.F), and the bridging scenario (III.G). Conclusions 
are summarized in Sec. IV. 

II. THEORY 

Following Muthukumar's work with monovalent salt^, we consider a linear flexible poly- 
electrolyte chain of N monomers in a solution of volume Q, with the center of mass of 
the chain at the origin of the coordinate system. Each monomer is monovalently charged 
(negative) and of length /. There can be either or both monovalent and divalent salts in 
the solution (say, water), which being electroneutral at all times will have a maximum of N 
monovalent counterions in addition to the salt ions. We assume that the counterion from the 
monovalent salt (say, Na + from NaCl) is chemically identical to the counterion from the poly- 
mer (say, Na + from NaPSS). Similarly, the coions from both types of salts are of the same 
species (say, Cl~ from NaCl and BaCl 2 ). At any time, both monovalent and divalent counte- 
rions (say, Ba ++ from BaCl2 as divalent counterions) can condense on separate monomers. In 
addition, the Ba ++ -monomer ion pair is viewed as a positive monovalent ion, and the negative 
coions (Cl~) will condense on some of these pairs as counterions. Therefore, if Mi monova- 
lent counterions and M 2 divalent counterions get adsorbed on the chain (Mi + M 2 < N), 
and M 3 (negative) coions condense on Ba ++ -monomer ion pairs (M 3 < M 2 ), the effective 
(or average) degree of ionization of the entire chain is / = (1 — Mi — 2M 2 + M 3 )/N. R g is 
the radius of gyration of the chain. c s \ and c s2 are, respectively, the number concentrations 
of the added monovalent and divalent salts. Both types of salts are fully dissociated into n\ 
monovalent counterions (Na + ), n 2 divalent counterions (Ba ++ ) and ni + 2n 2 negative coions 
(Cl~). Therefore, c s \ = ni + /Q and c s2 = n 2+ /Q. The free energy of the system, consisting 
of the chain, condensed and mobile counterions, and the solution, would depend on four 
independent variables: Mi, M 2 , M 3 , and R g . The theory^ aims to evaluate Mi, M 2 , M 3 , and 
R g self-consistently by calculating the free energy F of the system as a function of all these 
variables and electrostatic parameters. The equilibrium values of the four major variables 
are obtained by minimizing F simulatnously with respect to these variables. Therefore, this 
requires an extension of the previous theory of minimization with two variables to a theory 
with four variables. 

As before^, the free energy F has six contributions Fi, F 2 , F 3 , F 4 , F 5 , and F 6 related, 
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respectively, to (i) entropy of mobility of the condensed counterions and coions on the poly- 
mer backbone, (ii) translational entropy of uncondensed counterions and coions (including 
salt ions) which are mobile within the volume Q, (iii) electrostatic fluctuation interaction 
(Debye-Huckel) among all unadsorbed counterions and coions except the polymer, (iv) the 
unscreened electrostatic (Coulomb) energy of monomer-counterion pairs (both monovalent 
and divalent counterions) and monomer-counterion-coion triplets, (v) free energy of the poly- 
electrolyte with an average degree of ionization /, and (vi) electrostatic correlation involving 
the neutral ion-pairs or ion-triplets along the backbone of the polymer. 



A. The free energy 

To determine the entropic contribution from the condensed counterion and coions, we 
note that there are N monomers, Mi condensed monovalent counterions (Na + ), M 2 — M 3 
condensed divalent counterions (Ba ++ ) with no negative coion (CI - ') condensation, and M 3 
ion-triplets ('monomer-Ba ++ -Cl~'). Therefore, N — Mi — M 2 monomers remain bare charged 
with no ions condensed on them. Consequently, the partition function is 

N\ 

v — (o i) 

1 (N - Mi - M 2 )\Mi\(M 2 - M 3 )\M 3 \' y ' J 



We define, 



Mi M 2 M 3 

«i = ^; a2 = ^T' a3= AT (2 - 2) 



Using Fi = —KbTIiuZi, we have, 

Fi 

Nk rp = (! - a i ~ a 2) log(l - ai - a 2 ) + ailogai (2.3) 
+ (a 2 - a 3 ) log(a 2 - "3) + «3 loga 3 . 

The above expression implies two obvious constraints: 

Oil + a 2 < 1 an d a 3 < a 2 . (2.4) 

To determine the translational entropy of the uncondensed ions which are distributed in 
the bulk volume Q, we count as mobile ions: iV — Mi + rii + monovalent counterions (Na + ), 
n 2+ — M 2 divalent counterions (Ba ++ ), and ni + + 2n 2+ — M 3 monovalent negative coions 
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(Cl~). Therefore, the partition function related to the translational free energy in volume 
fl would be 

QN-M 1 +n 1+ +n 2+ -M2+n 1+ +2n 2+ -M3 



(2.5) 



(N — Mi + n 1+ )\(n 2+ - M 2 )!(n 1+ + 2n 2+ - M 3 )! 

QiV-Afi-Ai'2-Af3+2ni + +3n2+ 

~ (iV - Mi + n 1+ )\{n 2+ - M 2 )!(ni+ + 2n 2+ - M 3 )!' 
Relating numbers of ions and their concentrations is helpful. We note 

at n CslN Cs2N foa\ 
N = pU; m+ = ; n 2+ = . (2.6) 

P P 
With use of F 2 = —KbT In Z 2 and after some calculations we arrive at 

+ — ) log(p(l - ai) + c si ) + (— - a 2 ) \og{c s2 - pa 2 ) (2.7) 



Nk B T y p ' oyry ±J v p 

+ (— + 2— - a 3 ) log(c s i + 2c s2 - pa 3 ) 
P P 

- {(1 — ai — a 2 — a 3 ) + 2— + 3—}. 

P P 

Here, the constraint would be 

M 2 < n 2+ . (2.8) 

The free energy contribution from the correlations of all dissociated ions can be calculated 
in the electrostatic free energy 

F: ' = ~ kBT % (2 - 9) 

where the inverse Debye length k is given by 

A-kIb^Z^/VL. (2.10) 



K 2 



This result is obtained from the Debye-Huckel (DH) theory and caution must be exercised 
to identify regimes where the DH theory is at best a gross approximation. Here, Z± is the 
valency of the dissociated ion of the z-th species. In this case [see the text before Eq. (12. 5p ], 

k 2 = 4nl B {N - Mi + + 4(n 2+ - M 2 ) + n 1+ + 2n 2+ - M 3 }/Q (2.11) 
= AttI b (N - Mi - 4M 2 - M 3 + 2n 1+ + Qn 2+ ) /fi 
= 4nl B {p(l -a x - Aa 2 - a 3 ) + 2c sl + 6c s2 } . 
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Using the definitions from Eq. ( 12. 6j) we reach 

= _iV4^4 /2 -{p(l - ai - 4a 2 - a 3 ) + 2c sl + 6c s2 } 3/2 . (2.12) 

To determine the electrostatic energy gain due to condensation of all sorts of ions (Na + , 
Ba ++ , and CT~), we recount different ion pairs and triplets which form after condensation. 
On the polymer chain, there are ot\N pairs of 'monomer(-l) and Na + ion', (a 2 — 013) N pairs 
of 'monomer(-l) and Ba ++ ion', and a^N triplets of 'monomer(-l), Ba ++ and CI - ions'. In 
addition, (1 — a± — a 2 )iV monomers(-l) remain charge uncompensated. In Muthukumar's 
theory^, the dielectric mismatch parameter 5 was conceived to address a local dielectric 
constant e; in the vicinity of the chain backbone. Experiments have shown 33 q to be almost 
an order of magnitude less than its bulk value e (around 78 in water) near polyelectrolyte 
or protein backbone. The dielectric constant increases exponentially 34 - 35 from the fractional 
to its full bulk value over a distance of 1-loA from the chain monomers. 5 = (el/eid) was 
introduced^, where d is the dipole length of the monomer-monovalent counterion ion-pair. A 
glance at Fig. [2] reveals that 8 in above form only applies to the monomer-monovalent (Na + ) 
and monomer-divalent (Ba ++ ) ion pairs, but not to the divalent counterion- monovalent coion 
(Ba ++ -Cl _ ) ion pair in the monomer-divalent counterion-monovalent coion triplet. In an 
ion pair, there are two ions involved with a fixed distance between them. For the triplet, 
however, there are three lengths involved (for example, Ba ++ -monomer, Ba ++ -Cl _ and CT~- 
monomer), and interpretation of 8 is a bit tricky. We introduce a parameter 82 for the 
'monomer-Ba ++ -Cl - ' triplet. 5 2 is expected to be less than 48 (the value it would have 
assumed if there were two point charges, +2e and — 2e, respectively), but the determination 
of its actual value would probably require a microscopic treatment. In principle, 82 would 
be a function of 5. For simplicity, we assume all ions and monomers to be of the same size, 
and determine 8 2 in what follows. 

First of all, we write the electrostatic energy of condensation in terms of S and 8 2 '- 

F 4 ~ 
— — — = -atdls - 2(« 2 - a 3 )8l B - a 3 5 2 lB, ( 2 -13) 
JM kb± 

where Ib = l B /l and the terms containing 5 are written following Ref. 1 (with I ~ d). Note 
that if one assumes that the local dielectric constant t\ applies only to the Ba ++ -monomer 
pair, but not to the Ba ++ -Cl _ pair (e = e wa ter in that case), then <5 2 turns out to be 

62= (2 + -) 6. (2.14) 
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Those should be the lowest values for 82- On the other hand, if ei applies to both Ba ++ - 
monomer and Ba ++ -Cl~ pairs, then 

6 2 =46. (2.15) 

Those should be the highest values for 62- In practice, 82 would be somewhere between these 
two limiting sets of values. We choose the dielectric constant to be ei for the Ba ++ -monomer 
pair and (q + €b u ik)/2 for the Ba ++ -Cl~ pair. Then, 82 turns out to be 




(2.16) 



Point worthy of note is that the repulsion between the monomer and the Cl~ ion has been 
ignored; it would bring a very small correction in all three cases above. Although the 
counterion distribution and chain conformations are sensitively dependent on 62 , we later 
show that the very basic qualitative results do not change if 82 is assigned any value in the 
range mentioned above. Therefore, we would use Eq. (I2.16P in all of our representative 
calculations, unless mentioned otherwise. 

The free energy of the flexible polyelectrolyte chain is obtained by the variational 
method*^ in which one starts from the Edwards Hamiltonian, 

h = i/ Lrf K^) 2+ ^/ Lrfs r^ ( ^ (s) "^ (s/)) (2,i7) 

+ l fl ds l ds '\m^wy\ exp ~ Klnis) ~ nis% 

where L = Nl, lZ(s) is the position vector of the chain at arc length s, and w is the 
strenth parameter for all short-ranged hydrophobic or excluded volume effects. An effective 
expansion factor l\ is defined as follows: 

(R 2 ) = Nlh = Nl 2 h = QR 2 . (2.18) 

Here, (R 2 ) is the mean square end-to-end distance, and l\ effectively measures the swelling 
of the chain compared to a Gaussian chain. Assuming uniform, spherically symmertic ex- 
pansion or contraction of the chain, and by extremizing the free energy we obtain 
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where 
where 



)exp(a)e rf c(^) + i- + l-^|-^|, (2.20) 



a = K 2 Nlx/§. (2.21) 

Here, k = k/. We further define two more dimensionless variables, p = pi 3 and Csj = c s ,/ 3 , 
where i stands for the ion species. The important factor / is our previously defined average 
degree of ionization and is given by 

/ = 1 - ai - 2a 2 + a 3 . (2.22) 

The justification of using the variational result as the polymer free energy has been discussed 
in the previous paper-. We assume that similar arguments are still valid for this work. Any 
other alternative function of average degree of ionization (/) and radius of gyration (R g ) for 
F 5 may be used in place of Eq. (12.191) . 

Till now we have considered the electrostatic interaction only between monomers the 
effective charge of which, with or without condensed ions, is non-zero. In other words, the 
third-term in the polymer free energy F 5 [Eq. d 2 . 1 9 [) ] addresses the electrostatic interaction 
(within Debye-Hiickel approximation) between monomers with non-zero effective monopole 
charges. Further, F 5 considers only the monopole contribution of each ion-pair or ion- 
triplet. For example, a monomer-Na + pair and a monomer-Ba ++ -Cl~ triplet would be 
treated equally by F 5 , although they have quite different electrostatic effects. Similarly, a 
monomer-Ba ++ pair would be simply treated as a +1 charge, although the pair will have 
additional dipole effects. These additional dipole or higher order multipole effects would be 
critical when the average charge of the chain is close to zero. It has been shown^^! that 
these ion-pair effects play a key role to collapse a chain in presence of monovalent counterions 
at very low temperatures (i.e, when the degree of ionization is negligible). In the previous 
paper 1 , this correlation among neutral ion-pairs and bewteen neutral ion-pairs and charged 
monomers were addressed by short-ranged 5-function potentials which led to free energy 
contribution of the form 

4 ( 3 -) 3 "«>i6 2 '~B 2 ^^», (2-23) 
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where Wi's (< 0) are temperature dependent parameters, and are different for dipole-dipole 
and dipole-monopole interactions. These contributions are attractive and would modify the 
excluded volume interaction [the second term in F 5 , Eq. (12.191) ]. They can significantly 
reduce the size of the chain only around the isoelectric point (/ ~ 0), and the type of 
collapse is generally continuous or second order. 

In addition to the short-ranged dipole correlations, there can be long-ranged attraction 
between monomers mediated by multivalent counterions^ 1 ^. This attractive correlation 
between counterions may compensate the residual Coulomb repulsion of the chain even at 
higher degrees of ionization 40 , and the extended conformation of the chain may become 
unstable. This can as well be treated with the concept of ion 'bridging— It is still not 
conclusively known what kind of collapse this correlation-induced long-ranged attraction 
may induce. We leave out the short-ranged correlation effects near the isoelectric point in 
our present analysis. In Section III.G, we give preliminary results of an ion bridging theory 
leading to global instability of a polyelectrolyte chain based on our model. The point worthy 
of note here is that the bridging interaction reduces^ the effective value of the excluded 
volume parameter w. Therefore, for higher values of w, only very high Coulomb strength or 
divalent salt concentration will allow the bridging effect to take place. In most of our analysis 
(only except Section III.G), we assume w to be high enough to render the bridging effect 
to be negligible. Although we assume w to be zero except in Section III.G, that will imply 
the 'no-bridging' scenario in which choosing non-zero positive value of w only brings minor 
quantitative changes to our results. When bridging is included (Section III.G), however, w is 
a very important parameter affecting the transition salt concentration or Coulomb strength. 
A more detailed analysis of the role of multivalent cations in collapsing a polyelectrolyte 
and the related order of the transition will be presented in a future publication. 



III. RESULTS AND DISCUSSION 



We can express the total free energy F = F\ + F 2 + F3 + F4 + F5 in terms of the fraction 
of condensed counterions and coions (en's), size of the polymer (Zi), temperature and bulk 
dielectric constant (Ib), degree of polymerization (A), monomer density (p), monovalent and 
divalent salt concentrations (c s 's), and local dielectric mismatch parameters (5 and 5 2 ). The 
goal is to self-consistently determine the fractions of condensed ions (ai, a 2 , and a^) and the 
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size (R g = y/(Nlh/6)) that minimize the free energy. It is a simultanous minimization with 
respect to four variables («i, a 2 , a 3 , /i) instead of two in the previous paper-, and it is best 
performed numerically. Compared to a neutral system, there are two additional length scales 
in a charged system. They are the Bjerrum length (Ib) related to the Coulomb interaction 
and the Debye length introduced by screening due to dissociated ions including salt 

ions. This formalism invokes a third length scale due to the dielectric mismatch parameter 
8. Therefore, the important parameters on which we base our analysis are Ib, the salt 
concentrations (c a i and c S 2), and 8. 

A. Competitive adsorption 

We start with an isolated polyelectrolyte at low concentrations and at modest tempera- 
tures. In the first observation, the concentration of the divalent salt (c S 2) is increased from 
while keeping the concentration of the monovalent salt at zero (c s i = 0). We have chosen 
a higher and a lower value of 8 = 2.5 and 1.5, respectively. Generally at higher values 
of 8, fractions of condensed ion species are expected to increase. We notice that [Fig. [3] 
both divalent counterions and negative monovalent coions condense progressively in higher 
numbers (a.2 and a%, respectively) with increasing divalent salt concentration. The number 
of monovalent counterions («i), however, decreases with increasing c S 2- This implies that 
in this competitive adsorption process, condensed monovalent counterions, when challenged 
by a divalent salt, are replaced by divalent counterions. This happens for the entire physical 
range of the dielectric mismatch paremeter (we will later show in the diagrams of charged 
states). The variable values chosen in this specific calculation are degree of polymerization 
iV = 1000 and monomer density p = pi 3 = 0.0005 at Is = 3.0 (value related to flexible poly- 
mers of the type sodium polystyrene sulphonate (NAPSS) in water at room temperature). 
82 is given by Eq. ( 12. 16ft throughout the paper, unless noted otherwise. For 8 = 1.5, only 5% 
of monomers are neutralized by monovalnt counterions at no salt situation (c s i, c S 2 = 0) [Fig. 
E(a)] whereas the number increases to 35% for 8 = 2.5 [Fig. EJ^b)]. At this higher 8 value, al- 
most all available divalent counterions condense displacing the monovalent counterions with 
increasing c S 2- Negative monovalent coions (Cl~) also condense on the monomer-divalent 
ion-pair substantially. / decreases monotonically and reverses sign when c S 2 is still well 
below the monomer concentration p. At c S 2 ~ 80% of p, a± drops below 5% and «2 is about 
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80% implying that almost all available divalent counterions have condensed. The size of 
the chain decreases steeply [Fig. E](d)] due to rapid neutralization, but increases beyond 
the isoelectric point (/ = 0) due to repulsion among divalent cations that overcharge the 
chain. The size is dictated by the third term in Eq. (12.1 9p at these salt concentrations. For 
higher values of c S 2, more negative coions condense on the chain to reduce the (over) charge 
of the chain marginally. Number of condensed monovalent ions, however, decrease to zero 
monotonically. 

For the lower value of 5, the original sign of the polyelectrolyte charge (/) is preserved 
even at higher divalent salts [Fig. [3(a)] with the minimum absolute degree of ionization 
being around 0.27. Consequently, the size of the chain [Fig. [3(c)] remains substantially 
bigger than the Gaussian value for the entire range of salt cincentration. 

B. The chain collapse 

In Fig. E](b), we have noticed that the polyelectrolyte net charge due to condensation of 
divalent counterions becomes negligible as soon as the salt concentration reaches half the 
polymer concentration (c S 2 ~ p/2). Consequently, the chain collapses to its Gaussian size at 
around this isoelectric point. To compare with the case of monovalent counterions, we plot 
both the degree of ionization / and the expansion factor li at Ib — 3.0 in Fig. [H The other 
parameters are: N = 1000, 5 = 2.5, and p = 0.0005. We notice that for the monovalent 
salt, degree of ionization / of the polyelectrolyte decreases moderately and monotonically, 
and never changes sign. Consequently, the size (Zi or R g ) also decreases monotonically with 
the Gaussian statistics being obtained only at very high salt concentrations (or at very low 
temperatures). For the divalent salt, however, the isoelectric point is achieved as soon as 
there are sufficient number of divalent counterions available to neutralize the chain. That 
happens at a very low c s2 . As a result, the polyelectrolyte collapses (Gaussian statistics) 
near this isoelectric point. 

This collapse of a generic polyelectrolyte (NAPSS) in water occurs for modest values of 5, 
at a modest presence of the divalent salt, and at room temperature. This phenomenon has 
long been noticed theoretically^, and very recently in experiments^ 1 ^ 1 ^ and simulations^. 
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C. The issue of overcharging 

Our theory predicts that charge neutralization and subsequent charge reversal would 
occur to an isolated flexible polyelectrolyte in aqueous solutions at room temperature and at 
a modest presence of a divalent salt. The parameter 5 in our theory plays an important role in 
the charge reversal induced by counterion condensation. Temperature is also an important 
factor regulating the relative weight of electrostatic energy gain of ion condensation. To 
show these effects, we plot the degree of ionization / and the expansion factor l\ of the 
chain as functions of Bjerrum length Z_g (inverse temperature and bulk dielectric constant) 
for various 5 values in Fig. [5J The other parameters are: N = 1000, p = 0.0005, c a \ = and 
c~ S 2 = 0.0005 = p. The concentrations of the divalent salt and the polymer is chosen to be 
equal to ensure the availability of enough divalent ions to condense over every monomer if 
physical conditions permit. In Fig. 13(a), we notice that there is negligible condensation for 
5 = 1 (which is the comparable value of 5 in simulations^ 1 ^) . Similar to the monovalent case 
[Fig. 2(a) and 3 in Ref.-], the chain is neutralized only at very low temperatures (there is a 
factor of two in Is because the Coulomb energy gain for each ion-pair is twofold for divalent 
ions). At no temperature there is overcharging for 5=1. There is, however, a drastic 
qualitative change in the dependencies of / and R g on lg for 5 values 2 and above. At a 
particular temperature To, the chain is neutralized and if T is further reduced, overcharging 
occurs (and the chain swells). T is higher for higher values of 5 as expected (T" 1 and 5, 
both favor higher degree of condensation). The absolute value of maximum overcharge and 
re-swelled size increase with 5 as well. In particular, the reswelled size is larger than the 
original swelling for 5 = 4.0. This is despite the absolute effective charge being lower at 
the maximal reswelling because, at this point, the Coulomb strength for this large 5 value 
is high enough to have repulsion between monomers stronger than at point of maximal 
swelling as we increase Is- Another point of note is that for higher temperatures, just as 
for monovalent counterions, only a fraction of available divalent ions condense. The optimal 
temperature at which the chain reexpansion is maximum shifts to a higher value with higher 
values of 5. It might be instructive to note that we increase c S 2 in Fig. [31(b), by fixing the 
system at the abscissa value of 3 in Fig. [5] and staying on the 5 = 2.5 curve. For very low 
temperatures, sufficient number of negative coions (CI - ) condense to gradually re- neutralize 
the chain for all 5 values. 
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To further explore the issue of overcharging, we plot / and l\ against S for Ib = 3 in 
Fig. El The other parameters are the same as in Fig. [5j At room temperature in aqueous 
solutions, there can be no overcharging unless 5 > 1.7. Only for S as high as 1.7, the dieletric 
heterogeneity would be strong enough to electrostatically guide enough divalent ions that 
condense and reverse the charge of the chain. The strong sensitivity of the total charge and 
conformation of the polymer on 6 is manifest in Fig. El in which / decreases from about 
93% to zero (and subsequently l\ decreases from about 25 to 1 - the Gaussian value) for 
S changing only from 1-1.7. For very high values of 5, Cl~ ions condense progressively at 
higher numbers to reduce overcharging. 



D. Condensation of CI ions 



In a monomer-Ba ++ -Cl~ ion-triplet, the colinear arrangement of the three charges, in the 
same order as written here, is electrostatically the most favorable one. This is true regardless 
of the strength of S or the range of the local dielectric constant e/. We assume that all triplets 
in the system have this specific colinear arrangement in which the line joining the charges is 
perpendicular to the chain backbone [see Fig. [2]. We have discussed before the ambiguity 
in determining the electrostatic energy gain per triplet formation. The highest and lowest 
values permitted by physical conditions of the strength paramater 82 related to the formation 
of the triplet has been determined in terms of 5 through Eqs. fl2.14p - fl2.15p . Although 
suggested by a few authors^, we did not consider a different local dielectric constant for 
isolated Ba ++ or Cl~ ions. The dielectric constant relevant for the Coulomb interaction 
between the Ba ++ and Cl~ ions in Eqn. (12.161) is different than e. This is because of the fact 
that these two ions are in the local environment of the chain backbone, and consequently, the 
local dielectic behavior of the polyelectrolyte (i.e., 5) exclusively determines the value of 62- 
In this theory, to minimize the number of adjustable parameters, we assume that the sizes 
of the ions (Ba ++ , Cl~, and Na + ) are of the order of the size I of the monomer. To illustrate 
that these approximations do not compromise the generality of the problem, we plot the 
fraction of condensed ions (ai, «2, «3), degree of ionization (/), and the size expansion factor 
(/1) of the polymer as functions of 62 ■ We vary 82 from its lowest [Eq. (12.141) ] to its highest 
value [Eq. (12.151) ] for a specific value of 5. In Fig. El^a), we choose 5 = 2.5, for which 
overcharging is evident for Ib = 3.0. Other parameters are: N = 1000, p = 0.0005 = c~ S 2 
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and c~ s i = 0. We notice that with increasing 82 , progressively larger number of Cl~ ions 
condense on the chain. This reduces overcharging and close to the highest value of 82, 
overcharging would be marginally eliminated. This result has close resemblance to recent 
simulations^, where it is observed that smaller ion sizes (which leads to effectively higher 
82 values) reduce the degree of overcharging. In Fig. W(Jo), we choose 8 = 1.5, for which 
there is no overcharging at Ib = 3.0. We notice that a change in 82 has negligible effect on 
both / and R g in this case. However, fractions of both condensed divalent counterions and 
monovalent coions (a^ and 0J3, respectively) increase with 82 , keeping the overall degree of 
ionization (/) approximately unaltered. 



E. The free energy profile 

One of the advantages of our equilibrium adsorption theory is that it is possible to 
compare the contributions of different factors in the total free energy (i*\ to F 5 ) as functions 
of the critical parameters. The major conclusion of the theory^ has been that the equilibrium 
distribution of counterions and the size of the polyelectrolyte are determined essentially by 
the competition between the translational entropy of dissociated ions and the Coulomb 
energy gain of condensed ions. This is indeed borne out by our calculation in the presence 
of divalent salts too, as shown in Fig. [HJ In Fig. [H^a), the separate parts of the free energy are 
plotted against the Bjerrum length Ib for a fixed divalent salt concentration (c~ S 2 = 0.0005, 
equal to the monomer density) and for a specific strength of dielectric mismatch (8 = 2.5). 
The major contributions to the total free energy come from the translational entropy F2 
[Eq. (12.71) ] and the Coulomb free energy F4 [Eq. (I2.13p ], For higher temperatures (lower 
/_b's), the entropic term is favored as electrostatics remains negligible compared to thermal 
fluctuations. For lower temperatures, electrostatics becomes progressively relevant, and 
many ions condense reflecting substantial gains in F±. The entropic contribution F\ [Eq. 
( 12.31) ] related to the mobility of condensed ions along the backbone has negligible effect, and 
so does the Debye-Hiickel contribution F3 [Eq. (I2.12p ] at these salt concentrations. 

In Fig. E^b), similar free energy components are plotted against 8 at the same salt con- 
centration and for Ib = 3.0. The curves in (a) and (b) are remarkably similar demonstrating 
the equivalence of the parameters Ib and 8. According to this adsorption theory, reduction 
of any of temperature, the bulk dielectric constant e or the local dielectric constant e\ (near 
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the hydrophobic regions of the chain backbone) by a similar factor would induce very similar 
effects to polyelectrolyte behaviors. This is especially valid for modest values of Ib and S. 

F. The diagrams of charged states 

We have proposed earlier a tentative state diagram (Fig. [Tj) of the total charge / ( or the 
degree of ionization = 1 — a% — 2a 2 + 03) of the polyelectrolyte. In this subsection, we present 
the actual state diagrams calculated from our theory, as functions of three major variables 

- the Bjerrum length (Ib), divalent salt concentration (c S 2) and the dielectric mismatch 
parameter (5). In what follows, one of these variables is fixed and the diagram of states 
(regions of negative and positive degree of ionization) is calculated numerically as functions 
of the other two. Figs. l9]fT0l describe the complete limiting charged states, parts of which 
have already been discussed in detail in preceding subsections. 

In Fig. 0(a), the calculated state diagram at Ib — 3.0 is presented as a function of the 
divalent salt concentration c~ S 2 and the dielectric mismatch parameter 5. The state diagram 
is qualitatively similar to the proposed one [Fig. [TJ, with the strength of the Coulomb 
interaction being represented by 5 [Eq. (12. 13ft ]. and with Ib being fixed. To explain the 
diagram, we first choose a specific value of 5 = 2.5 (see Figs. EJ^b) andH]) and monitor the 
charged state with increasing divalent salt concentration. For low salt, there are not enough 
divalent counterions (Ba ++ ) to neutralize the chain and the polyelectrolyte preserves its sign 
of charge (state A) of salt-free conditions. At around c S 2 ~ p/2, which is half the monomer 
concentration, the charge of the polymer becomes zero (on the locus of first isoelectric points 

- the solid line). If c S 2 is increased further, the polymer charge is reversed (state B), and at 
around c S 2 ~ p, almost all monomers are neutralized by divalent counterions. The charge 
reversal is maximum at around this point (on the locus of maximum overcharging points - 
the dotted line). With c S 2 increasing even further, more negative coions (CI - ) are available 
in the solution and some of them condense on the monomer-Ba ++ ion-pairs to reduce the 
degree of overcharging (state C). The first isoelectric points between states A and B are 
reached at a higher c S 2 for a lower 5, because a higher fraction of divalent counterions would 
remain dissociated in the solution due to a lower Coulomb energy gain. For values of 5 
higher than ~ 3, a substantial fraction of monovalent counterions of the polymer (Na + ) 
too remain condensed on the chain and the chain charge is neutralized with fewer divalent 
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counterions. If 5 is less than ~ 1.7, the state of overcharging (state B) is never reached 
and with increasing salt concentration, the polymer charge goes through a minimum (on 
the locus of points of minimum charge - the dashed line) before increasing again due to CI - 
ion condensation. The line of minimum charge (for S less than ~ 1.7) expectedly continues 
to be the line of maximum charge reversal (for 5 greater than ~ 1.7). For very high salt 
concentrations, the Coulomb interaction is progressively screened and all condensed ions 
begin to rejoin the solution (not included in the state diagram). We must, however, be 
cautioned that the Debye-Hiickel (DH) approximation (and consequently _F 3 in Eq. (12.121) ) 
might not be valid at this high salt regime. A salt concentration for which the Debye length 
is equal to the Bjerrum length (k _1 > Ib) can be tentatively set as the highest limit of 
validity of the DH theory. For divalent salts it turns out to be [see Eq. (12.111) 

C s2 -max ^ (247r/|) _1 . (3.1) 

The steepness of the state boundary (the locus of second isoelectric points - the dot-dashed 
line) implies that the polymer charge becomes zero again (only applicable for 5 > 1.7) due 
to re-dissolution of condensed ions at least an order higher salt concentrations. 

The state diagram as a function of Ib and c S 2 for a fixed value of 5 = 2.5 is presented in Fig. 
E^b). The diagram is qualitatively similar to the previous one, although a much higher salt 
concentration (note the difference in the scale of the coordinate) is needed to reach the line 
of minimum charge (dashed) and the line of maximum overcharging (dotted) at low values of 
Ib (higher temperatures). In this regime, electrostatics becomes progressively weaker with 
increasing temperature, and consequently lower fractions of available divalent ions condense. 
Regarding this diagram too, the degree of ionization and overcharging (absolute value of /) 
can be obtained for the particular value of Ib = 3.0 from Fig. E](b). 

The state diagram as a function of Ib and S for a fixed salt concentration c S 2 = 0.0005 
(equal to the monomer concentration) is presented in Fig. [101 One essential characteristic is 
that the degree of ionization remains steadily at zero (state CI) above a certain value of Ib 
(i.e., below a certain temperature), because the Coulomb attraction is strong enough to make 
form the monomer-Ba ++ -Cl" ion-triplet on every monomer location. This critical value of 
Ib (on the lower boundary of the zero charge state - the dot-dot-dashed line) decreases 
with higher values of 5 (higher electrostatic energy gain). The magnitude of the degree 
of ionization as a function of Ib for fixed values of 5 can be obtained in Fig. [51 and as a 
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function of 5 for a fixed value of Ib in Fig. [6j Both figures can be analyzed in conjunction 
with this state diagram. One small point worthy of note is that state D, although implying 
a non-zero degree of ionization of the same sign of the bare polymer charge, has virtually 
negligible charge in it that can be assumed to be zero. 



G. The bridging scenario : a simple theory 

In all previous discussion, we did not consider the bridging configuration of nonbonded 
monovalent monomers by divalent counterions. This ion-bridging phenomenon is electro- 
static in nature and can significantly affect polyelectrolyte conformation if present. In this 
subsection we will give preliminary results of a simple theory based on our model. We as- 
sume that a fraction of condensed divalent ions participate in bridging. Therefore, when 
bridging is included the minimization of the free energy, as it will turn out, is with respect 
to five variables. Now, the bridge formed by divalent counterion tantamounts to a cross- 
link junction of functionality four—, that in turn can be treated as an attractive two-body 
interaction of local nature (like two-body excluded volume interaction). Therefore, in the 
presence of bridging effects due to divalent counterions, w in Eq. (12.191) is replaced by, 

w' = w + -r^a 2b , (3.2) 
k B T 

where a 2 b is the ratio of the number of divalent ions that participate in bridging to the 
number of monomers (i. e., a 2 b = M 2 b/N, where M 2 b is the total number of divalent ions 
involved in bridging). E br is the attractive energy associated with one bridge, and hence is 
negative. To calculate Eb T , the relevant dielectric constant should be the local one (q) since 
the divalent cation in the monomer-cation-monomer charge complex sits between and in the 
vicinity of both monomers. With the definition 1 5 = el/ed, where the distance between both 
ion pairs in the complex remains to be d — I (same as other condensed pairs), it turns out 
that 

4e 2 e 2 
Ebr = —~ ; + 



= - 7 -l B 5k B T. (3.3) 

We must add the third virial term in the chain free energy (F5, Eq. (12. 19ft ) to maintain 
stability in the system in case of a negative w'. Combining Eqs. (13. 2\\ and (13.31) . therefore, 
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the chain free energy takes the form, 



4 ~ 1 " logO + t f^:V /2 (™ " h~BSa 2b )4=^L (3.4) 



Nk B T 2N \ * L J 3\2nJ v 2^ 'v^V^ 3 / 2 

lmj in /6 , 2; ~ iV 1/2 , 



where W3 is the third-virial coefficient which is necessarily positive. Further we note that 
if a fraction «2b of condensed divalent counterions participate in bridging, a fraction a 2a = 
&2 — oi2b does not. Therefore, the electrostatic energy related to the formation of monomer- 
cation monocomplexes (F4, Eq. ( 12.131) ) is modified after the inclusion of bridging interaction 

as, 

F A ~ 
— — — = -a x 5l B - 2(a 2a - a 3 )6l B - a 3 6 2 l B , (3.5) 
JM k B L 

where a 3 < a 2a . The other parts of the free energy remain unaltered and are given as: F\ 
in Eq. ( 12. 3D . F 2 in Eq. ( 12.71) . and F 3 in Eq. ( 12.121) . In all these cases, a 2 = a 2a + a 2 b- 

The total free energy F = F% + .. + F5 is minimized now for a new set of five variables, 
«i, a 2a , oc 2 b, «3 and li, and the polymer and counterions are free to explore every possible 
degree of freedom. The representative result is given in Fig. [TTJ The parameters chosen 
are: N = 100, p = 0.0008,4 = 3.0,5 = 1.9, c sl = 0,w = 2.0, w 3 = 0.25. For very low 
divalent salt concentrations, the conformations are very similar to the case where bridging 
is absent. At modest temperatures (l B — 3.0 in water) and for low salt (c~ s2 < 0.00027), 
almost all added divalent counterions condense, but they from monocomplexes (no bridging, 
Fig. [TTTa)). At a particular c~ s2 *, which depends on the prevalent physical conditions, all 
divalent ions suddenly form dicomplexes (bridging) accompanied by a collapse of the chain 
[Fig. fTTTb). in which li <C 1 for c~ s2 > 0.00027] and a huge gain in electrostatic bridging free 
energy [Fig. [TTT c)]. The Cl~ ions condense as they do for no-bridging scenario only if the 
divalent salt concentration is lower than the collapse concentration c~ s2 * . Above that, Cl~ 
ions become free as every condensed divalent cation is attached to two monomers. The effect 
of the excluded volume parameter w is evident in Eqs. (13. 2p and (13.4)) as they show that 
a higher w will require a higher Coulomb strength or divalent salt concentration to effect 
the bridging collapse. Until the collapse, the distribution of counterions and the polymer 
conformations are quite similar to that of the 'no-bridging' cases (see Fig. [3] for example). 
This explains our choice of w = for the rest of the article (except for this subsection). In 
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the 'no-bridging' scenario, different positive values of w would only render minor qualitative 
changes to our results. 

We further notice that the collapse concentration, c~ S 2*, decreases with increasing Coulomb 
strength [Fig. [T27 a).(b)] confirming that the first order collapse induced due to ion-bridging 
by divalent (or multivalent) cations is an electrostatic phenomenon. In addition, Fig. [T2Tc) 
shows that c~ S 2* roughly varies inversely with both forms of Coulomb strength, Ib and 5. This 
is a remarkable prediction for experimentalists and we find Ibc~ S 2* — 0.0006 (for 5 = 2.5) and 
Sc~ s2 * ^ 0.0005 (for f B = 3.0). 

In conclusion, our model predicts a bridging transition, which we believe depends sensi- 
tively on temperature and dielectric heterogeneity, as well as on the availability of divalent 
counterions. At this point, we leave the rest of the bridging analysis for future publication. 

IV. CONCLUSIONS 

We have extended Muthukumar's adsorption theory 1 for condensation of monovalent 
counterions on a flexible polyelectrolyte by including divalent counterions. It is observed 
that the divalent counterions replace the monovalent ones in the competitive adsorption 
process. For moderate values of the dielectric mismatch parameter (6), temperature and bulk 
dielectric constant (Ib), numbers of both condensed divalent cations (Ba ++ ) and monovalent 
anions (Cl~) increase and the number of condensed monovalent cations (Na + ) decrease 
monotonically with increasing concentration of added divalent salt (BaC^). As observed in 
previous theories, experiments, and simulations, a moderate amount of divalent salt entirely 
compensates the polyelectrolyte (NaPSS) charge and consequently contracts the chain to its 
Gaussian size. The divalent salt concentration at which the charge neutralization takes place 
is roughly half the salt concentration implying that almost all divalent cations added to the 
solution condense at those modest parameter values. With further increase of divalent salt 
concentration, condensed divalent counterions overcharge the chain resulting in reswelling of 
the chain. This phenomenon is addressed for the first time theoretically for a flexible chain 
allowed to take all possible conformations. The charge reversal, however, is absent at modest 
temperatures regardless of the salt concentration if the dielectric heterogeneity is weaker than 
a critical value. The minimum degree of heterogeneity (5) required to condense all available 
divalent cations at a particular temperature and solution (Ib) increases with temperature. 
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At room temperature in aqueous solutions (Is ~ 7A), the adsorption theory predicts that 
the ratio of the bulk to local dielectric constant must be at least 1.7 to overcharge the 
polyelectrolyte. 

It can be conceived that there are three length scales involved in the polyelectrolyte 
system analyzed by this adsorption theory. They are, the Bjerrum length lg (representing 
the equilibrium temperature T of a solution with bulk dielectric constant e) , the Debye length 
hT x (representing the salt concentration c s i), and the third one related to the strength of 
the dielectric heterogeneity in the vicinity of the polyelectrolyte backbone (5). We have 
developed the diagrams of charged states of the polymer in terms of these three variables. 
Both higher 6 (lower ratio of local to bulk e) and lg (lower temperature or bulk e) facilitate 
ion condensation of all types. Consequently, the isoelectric point (charge neutralization) and 
subsequent charge inversion are achieved at lower temperatures for higher 5 and vice versa. 
The absolute value of the maximum charge reversal and simultaneous maximal reswelling 
increase with 5. For a minimal dielectric mismatch (5 ~ 1), no overcharging is predicted 
for any temperature. This limit is the closest comparable to the simulations and it explains 
why in some simulations overcharging is never observed^. 

For moderate values of 8 (> 2), maximal overcharging and reswelling occurs at intermedi- 
ate temperatures lower than that of the first isoelectric point (at which a moderate amount 
of divalent cations and a fraction of monovalent anions neutralize the chain). If temperature 
is further decreased (Ib increased), progressively higher fraction of monovalent anions (CI - ) 
condense reducing the degree of overcharging. At low enough temperatures, regardless of 
the degree of dielectric heterogeneity, the second and trivial isoelectric point (at which all 
monomers are neutralized by the presence of both Ba ++ and Cl~ ions) is achieved and the 
Gaussian statistics is reestablished. 

We have determined the physical limits of the strength parameter for the formation of 
monomer-Ba ++ -Cl~ triplets as functions of the basic dielectric mismatch strength 5. It has 
been shown that the qualitative picture of the total effective charge and conformation of the 
polyelectrolyte remain unaltered although the distribution of counterions and coions around 
it may vary for a change in the triplet strength parameter (62) within these limits. For an 
unrealistic value of the parameter closer to its higher limit, increased Cl~ condensation may 
marginally eliminate overcharging at all physical conditions. 

Analysis of the free energy hightlights the competition between the entropic contribution 
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of dissociated ions and the electrostatic contribution of condensed ions in shaping the coun- 
terion distribution around the polyelectrolyte. The local dielectric constant and temperature 
have similar contribution to the electrostatics of the system (except at very high salt con- 
centrations). If electrostatic bridging of non-bonded monomers by divalent counterions is 
effective due to a weak excluded volume effect, there is a first order collapse at modest con- 
ditions. The transition salt concentration varies approximately inversely with the Coulomb 
strength at modest conditions. In summary our theory predicts that both electrostatics and 
the dielectric inhomogeneity near the backbone of the polymer (reflecting the nonuniversal 
chemical nature of the system) are responsible for charge reversal in polyelectrolyte systems. 
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FIGURE CAPTIONS 

Fig. 1 Sketch of charged states for an isolated polyelectrolyte chain (of NaPSS type) in 
dilute solutions (water) in the presence of a divalent salt (of BaCl 2 type) as functions 
of the Coulomb strength (l B ~ (eT) _1 ) and the salt concentration. Points left to the 
isoelectric line (on which net effective charge (degree of ionization) on the polymer is 
zero) correspond to states in which the sign of polymer charge is unchanged (negative). 
However, there is a locus of points for intermediate values of salt concentration at which 
the net charge is a minimum. Right to the isoelectric line the effective polymer charge 
is reversed (positive). If the isoelectric point is crossed along the line of minimum 
charge from left to right, it becomes the line of maximum overcharging. Dashed part 
of the isoelectric line is beyond the Debye-Huckel limit. 

Fig. 2 Schematic diagram of the system consisting of the isolated polymer chain, con- 
densed counterions, dissociated mobile ions, and the solution as the background in- 
teracting only through the dielectric constant e. Possible charge complexes for each 
monomer: monomer (-1), monomer-monovalent (-1,+1), monomer-divalent (-l,+2), 
and monomer-divalent-monocoion (-l,+2,-l). The dielectric constant e\ in the vicinity 
of the chain is much lower than the bulk value. To reach equilibrium the major com- 
petition is between the translational entropy of the dissociated ions and the Coulomb 
energy gain of the condensed ions. 

Fig. 3 Competitive displacement of monovalent counterions by divalent counterions: frac- 
tion of condensed ions (aii, a 2 , a 3 ), degree of ionization (/ = 1 — cei — 2a 2 + c* 3 ) 
in (a) and (b) for 5 = 1.5 and 2.5 respectively, and the expansion factor (/i = 



6Rg/Nl 2 ) in (c) and (d) for the same 5 values plotted against the divalent salt 
concent rat ion (c^2 ~ 0.01c s (M)). For Gaussian chain, Zi = 1. Other parameters are: 
N = 1000, p = 0.0005,4 = 3.0, w = 0, and c~ sl = 0. Note, for lower 5, there is 
no overcharging. For higher 5, almost all divalent counterions condense on the chain 
replacing the monovalent ones. The number of condensed negative coions (a^) closely 
follow «2 for this particular value of 5. The sign of / is reversed (overcharging) at 
some concentration of BaCl2. Near the isoelectric point (/ ~ 0), the chain is Gaussian 
due to minimal electrostatic repulsion. It swells due to overcharging if c~ S 2 is further 
increased. 

Fig. 4 Effect of valency of counterions: comparison of degree of ionization (/) in (a), and 
size expansion factor (li) in (b), of the polyelectrolyte in presence of either monovalent 
or divalent salt. 5 = 2.5 and other parameters are the same as in Fig. [3J c~ S 2 — 
is zero when c sl is added and vice versa. Divalent counterions can neutralize and 
consequently collapse the polymer at moderate conditions of Ib = 3.0. For monovalent 
counterions, collapse (to Gaussian chain) is only possible at very low temperatures. If 
c~ S 2 is increased beyond the isoelectric point, the chain expands due to overcharging. 

Fig. 5 Dependency of overcharging on Ib' degree of ionization (/) in (a), and size expansion 
factor (li) in (b), of the polyelectrolyte plotted against Ib for different values of 5. 
Parameters are: N = 1000, p = 0.0005, c s \ = 0, and c~ S 2 = 0.0005. Collapse and 
subsequent overcharging occur for higher values of 5. Isoelectric point is reached at 
lower Ib for higher values of 5. This overcharging behavior contrasts with Figs. 2(a) 
and 3 in Ref «• for monovalent salts. 

Fig. 6 Dependency of overcharging on 5: fraction of condensed ions («i, 02,03), degree 
of ionization (/) in (a), and the expansion factor (Z x ) in (b), plotted against the 
(^-parameter for Ib = 3.0. Other parameters are the same as in Fig. [5l Charge 
neutralization occurs for 6 about 1.7; overcharging is possible only if 6 > 1.7. CD 
condensation (013) is higher for higher values of 5, and that reduces overcharging. 

Fig. 7 Effect of dielectric mismatch for the coion: fraction of condensed CD ions (0:3), 
degree of ionization / in (a) and (b) for 5 = 2.5 and 1.5 respectively, and expansion 
factor (Zi) in (c) and (d) for the same S values plotted against 82 for its whole physical 



range possible for the respective value of 5. Other parameter values are the same as 
in Fig. [6j For higher 5 = 2.5, progressively higher values of 62 reduce and finally 
eliminate overcharging (and consequently the reexpansion of the chain) by increasing 
a 3 . For 5 = 1.5 (no overcharging at any temperature), fraction of both condensed 
divalent counterions and monovalent coions, a<i and a.3 respectively, increase to leave 
/ (and R g ) approximately unchanged. 

Fig. 8 Contributions to free energy: separate parts of the free energy for a fixed divalent salt 
concentration (c~ S 2 = 0.0005) as functions of (a) Bjerrum length Ib and (b) dielectric 
mismatch strength (5). For (a), 5 = 2.5 and (b), lg = 3.0. Parameters are: N = 
1000, p = 0.0005, c~si = 0. Energies are: Fi (dot) = entropy of mobility along the 
chain, F 2 (dash) = translational entropy of mobile ions, F 3 (dot-dash) = Debye-Hiickel 
correlation between mobile ions, _F 4 (dot-dash-dash) = Coulomb attraction between 
condensed ions and F tot (solid) = total free energy. For fixed salt concentration, the 
major competition is between the translational entropy (increase with temperature) 
and Coulomb attraction (increase with both Ib and S). The similarity in the roles of 
Is and S is evident. 

Fig. 9 (a) The state diagram of the total charge density of the polymer (/) for Ib = 3.0 
as functions of the dielectric mismatch 5 and the divalent salt concentration c~ S 2- Pa- 
rameters are: iV = 1000, p = 0.0005, c s \ = 0. Charged states are: A,D - negative, B,C 
- positive (note: original polymer charge is negative). Lines are: isoelectric branch 
one (solid), maximum overcharging (dot), isoelectric branch two (dot-dash), and min- 
imum charge (dash); (b) the state diagram of / at a fixed dielectric mismatch strength 
(5 = 2.5) as functions of the Bjerrum length Ib and the divalent salt concentration c~ S 2- 
Other parameters, states, and lines are the same as in (a). 

Fig. 10 The state diagram of the total charge density on the polymer (/) at a fixed divalent 
salt concentration (c~ S 2 = 0.0005) as functions of the dielectric mismatch strength 5 
and the Bjerrum length Ib- Other parameters are the same as in Fig. 9. Charged 
states are: A,D - negative, B,C - positive, and CI - zero. Lines are: isoelectric branch 
one (solid), maximum overcharging (dot), isoelectric branch two (dot-dash), minimum 
charge (dash), and zero charge (dash- dot- dot). 



Fig. 11 Effect of bridging: degree of ionization (/), fractions of condensed ions (ce's) (a), 
size expansion factor (li) (b), and separate parts of the free energy (F's) (c) as functions 
of divalent salt concentration (c S 2) when monomer-bridging by divalent cations is taken 
into account. Parameters are: N = 100, p = 0.0008, Ib = 3.0,5 = 1.9, c~ sl — 0, w — 
2.0, u>3 = 0.25. Bridging induces a first-order collapse transition with a sudden gain in 
electrostatic ion-bridging energy. At the transition salt concentration, all monomer- 
divalent cation ion-pairs ('monocomplex'es) give way to monomer-cation- monomer ion 
bridges ('dicomplex'es). 

Fig. 12 Effect of Coulomb strength on collapse: the salt concentration c~ S 2* \ at which the 
first-order collapse occurs, as functions of l B (a), and 5 (b). In (a), 5 = 2.5. In 
(b), Ib = 3.0. All other parameters are the same as in Fig. 11. Lowering of c S 2* 
with Coulomb strength indicates that the collapse is due to electrostatic interactions. 
/ b and 5 play similar roles, as expected. In (c): IbC S 2* — 0.0006 (for 5 = 2.5) and 
5c s2 * ^ 0.0005 (for l B = 3.0). 
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